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ABSTRACT 

We have conducted a differential elemental abundance analysis of unprecedented accuracy (~ 0.01 
dex) of the Sun relative to 11 solar twins from the Hipparcos catalogue and 10 solar analogs from 
planet searches. We find that the Sun shows a characteristic signature with a « 20% depletion of 
refractory elements relative to the volatile elements in comparison with the solar twins. The abundance 
differences correlate strongly with the condensation temperatures of the elements. This peculiarity 
also holds in comparisons with solar analogs known to have close-in giant planets while the majority 
of solar analogs found not to have such giant planets in radial velocity monitoring show the solar 
abundance pattern. We discuss various explanations for this peculiarity, including the possibility that 
the differences in abundance patterns are related to the formation of planetary systems like our own, 
in particular to the existence of terrestrial planets. 

Subject headings: Sun: abundances — solar system: formation — stars: abundances — planetary 
systems 



1. INTRODUCTION 

Studies of extra-solar planetary systems and the pos- 
sibility of life in the Universe depend funda mentally on 
whether the Sun and its planets are typical (iGustafsswil 
[200l iBeeTet all 120041: iRobles et al.l[200l ISeaeeri 12003) . 
Previous chemical composition studies have con c luded 
that the Sun is a typical star (iGustafssonj 1 19981 120081 : 
IRobles et al.l 120081 ). Besides an increased likelihood for 
stars with higher overall metallicity of hosting giant 
planets, no definitive chemical differences have been es- 
tablis hed between stars with and without known plan- 
ets (Gonza leslll997tlEcuvillon et al.l[2006l : lUdrv fc Santos! 
[1)07). Previous studies on chemical abundance anoma- 
lies in the Sun as compared with solar-like stars have 
been inconclusive due to the relativ ely large 0.05 
dex) remaining s ystematic err ors (iGustafssonl 120081 : 
IRobles et all [200l iReddv et al.] [20031). In such stud- 
ies, the selection of stars to which the Sun is compared 
is crucial. Different categories of solar-like stars can 
be identified. "Solar-type stars" range from late F to 
early K, "solar analogs" include only G0-G5 dwarfs, and 
stars almost identical to the Sun are called "solar twins" 
(|Cavrel de Strobel 119961 ). In modern stellar chemical 
analyses, the error budget is dominated by systematic 
errors in the model atmospheres and modeling of stellar 
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spectra (|Asplundll2005l ) . The highest possible accuracy in 
the analysis of the Sun relative to stars can be achieved 
in comparisons relative to solar twins, since the model 
errors will cancel to a large extent and effects of possi- 
ble systematic errors in temperature scale and absolute 
abundances are minimized. 

2. OBSERVATIONS AND ANALYSIS 

Our sample comprises 11 solar twins and 10 solar 
analogs. The solar twins were selected from more 
than 100,000 stars in the Hipparcos catalogue by em- 
ploying precise colour-temperature relations, trigonomet- 
ric parallaxes and age indic ators (jMelendez et al.l [20061 : 
iMelendez fc Ramirez I [20071 ). They lie within a distance 
of 75 pc from the Sun and belong to the Galactic thin 
disk population. No giant planets are known around our 
solar twins because most of them have not been searched 
for planets. The solar a nalogs were selected f rom radial 
velocity planet surveys (|Udry fc Santosl[2007f ). The Sun 
(reflected light from the asteroid Vesta), solar analogs 
and solar twins were observed with the MIKE spectrom- 
eter at the Clay 6.5 m Magellan telescope at Las Cam- 
panas Observatory in Chile, at high spectral resolution 
{R = AA/A — 65, 000) and very high signal-to-noise ra- 
tio {S/N ~ 450 per pixel at 600 nm) through the wave- 
length range 340-1000 nm. Example Magellan spectra 
of solar twins are presented in Fig. 1. Finally, one of 
our twins (HIP79672=18 Sco) was also observed with 
the HIRES spectrograph on the Keck I telescope to yield 
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an even higher quahty spectrum with S/N « 450 and 
R = 100, 000; a solar spectrum (asteroid Ceres) was also 

obtained with identical setup. 

A model independent ana lysis ([Melendez et al.l [20061 : 
IMelendez fc Ramirez I [2007*) indicates that the twins 
have effective temperatures within 75 K of the Sun, 
logarithmic surface gravities within 0.10 dex and loga- 
rithmic iron abundances [Fe/H] within 0.07 dex (where 
[A/B]= logioiNA/NB)star - logw{NA/NB)Q). These in- 
dications were confirmed by a detailed model-atmosphere 
analysis. The final adopted effective temperatures were 
obtained from the excitation equilibrium of Fe I, while 
the surface gravities were determined from the ionization 
equilibrium of Fe I/Fe II. We have achieved an accuracy 
in [X/Fe] at the level of 0.01 dex for the solar twins (e.g., 
the star-to-star scatter in [Cr/Fe] is only a — 0.009 dex). 
For the solar analogs high precision abundances are also 
obtained due to their similarity to the Sun, although not 
as high as for the solar twins. The observed scatter in 
relative abundances is fully consistent with the predicted 
error from the quality of the spectra and the remaining 
uncertainties in the stellar parameters. 

3. RESULTS 

Our results reveal that the solar chemical abundances 
relative to iron depart from the mean abundance ratios 
in the solar twins (Fig. 2). The Sun is enriched by ~ 0.05 
dex in most of the lighter elements, while other elements 
like aluminium are deficient by 0.03 dex, thus the off- 
set between the elements that are depleted and enhanced 
is ~ 0.08 dex (20%). In Fig. 3 the abundance differ- 
ences are s hown versus t he 50% condensation tempera- 
ture Tcond ([Loddersll2003l ) for a solar-system composition 
gas. Volatile elements (low Tcond) are more abundant in 
the Sun relative to the twins while elements that eas- 
ily form dust (elements with high Tcond, i-e. refracto- 
ries) are under-abundant. A very similar pattern, al- 
beit of much greater amplitude, is se en in the interstellar 
medium due to dust condensation ([Savage fc SembadJ 
[19.96), as well as in A Bootis stars and some highly 
evolved stars through th e accretion of dust-cleansed gas 
(IVenn fc LamberU [l99l . The correlation with Tcond 
shown in Fig. 3 is highly significant: the Spearman cor- 
relation coefficient is rs = —0.91 and an arrangement 
of this type have a probability of ^ 10~^ to happen by 
chance. The scatter around the mean trend in Fig. 3 
is only 0.01 dex. On a star-by-star basis, only 1-3 solar 
twins resemble the Sun chemically (Fig. 4). In fact, the 
probability of obtaining the same order between abun- 
dances differences for the four degrees of volatility for 9 
(and very nearly 10) out of 11 stars as a result of random 
errors is only 2x 10~^^. Previous attempts to trace Tcond- 
dependent abundance differences for stars known to host 
giant planets hav e not convincingly de monstrated signif- 
icant differences ([Udrv fc Santosll2007D . probably due to 
the high accuracy needed. 

The results for the solar analogs are shown in Fig. 5. 
The solar analogs with detected giant planets have an 
abundance pattern similar to the solar twins, while the 
solar analogs without detected giant planets are more 
similar to the Sun. 

In order to confirm that our results are robust, we have 
checked our abundance scale using spectra obtained with 
another telescope and spectrograph. We repeated our 



highly differential abundance analysis for the Sun and 
the solar twin HIP 79672, using a HIRES/Keck spec- 
trum. Our Keck-based results confirm the separation 
between volatile and refractory elements found with the 
MIKE data. The [X/Fe] abundance ratios of the HIRES 
analysis are in excellent agreement with our MIKE re- 
sults. In particular, the relative abundances of volatiles 
like carbon, oxygen and zinc, agree to within 0.003 dex 
between the HIRES and MIKE data. This demonstrates 
that our abundance scale is robust, and that the quality 
of the MIKE spectra is high enough to allow our conclu- 
sion that the refractories in the Sun are depleted with 
respect to the volatiles. 

We are not aware of any systematic errors that would, 
if corrected for, significantly reduce these differences be- 
tween the Sun and the twins. An error in the tempera- 
ture scale would lead to even larger departures in [X/Fe] 
but not generate the consistent pattern found here. 

4. THE ORIGIN OF THE SOLAR DEPARTURES 

In spite of the suggestions above that the solar charac- 
teristics could be related to dust condensation, we shall 
here first discuss other possible explanations. 
(a) Galactic- evolution effects: Perhaps the solar abun- 
dance pattern refiects the Galactic chemical evolution 
and the varying composition in the Milky Way ISM when 
the Sun was formed. The [a/Fe] ratio (a representing O, 
Ne, Mg, Si, Ca, S, Ti) was higher earlier in the history 
of the Galactic disk, due to the increasing role of SN la 
relative to that of SN II in enriching the ISM. The solar 
twins have iron abundances similar to solar. If the metal- 
licity of the Sun is unusually high for its age, most twins 
could be younger and have lower \a/ Fe\. The Sun is , 
howe ver, not so meta l-rich for its age (Fuhrmann '2004|; 
Hol mberg et al.l[2009f ) and it is unlikely that the variation 
of [a /Fe] with metallicity is large enou g,h for this hypoth- 
esis to work (iBensbv fc Feltzinid[200alChen et"aLll2002l : 
iNissen et al.ll2004D . Besides, the abundances relative to 
Fe for the volatile elements C, S and Zn as well as the re- 
fractory element Si all show no trend with metal content 
around solar metallicity. One direct way of testing this 
hypothesis would be to compare the ages of the twins to 
the Sun. Available niethods for estimating stellar ages 
([Melendez et al.l 120061: IMelendez fc Ramirez I I2007D ap- 
plied to the twins reveals that their median age is 4.1 
Gyr. Thus, the peculiar abundance pattern of the Sun 
is not the result of large age differences between the Sun 
a nd the twins. 

iWielen et al.l ([1996D suggested that the Sun may have 
migrated from an inner Galactic orbit. If there were ra- 
dial abundance gradients, systematic differences between 
the Sun and the twins may arise. However, the observed 
gradients in [0 /Fe] for youn g thin disk objects are not 
steep enough (iPrzvbillal 12008.) . and G alactic chemical 
evolution models ([Chiappini et al.l [2^0011 ) do not suggest 
any substantially steeper gradient 5 Gyr ago. 
(h) Supernova pollution: Another possibility is that the 
proto-solar cloud happened to be more polluted by a lo- 
cal supernova than other star formation regions in the 
Galaxy, although this hypothesis does not explain why 
the Sun, and not the twins, happens to be affected. How- 
ever, the peculiar solar abundance pattern cannot be 
explained by an offset of SNII relative to SNIa yields 
(,Wooslev fc'Weaver..l995; .Thielemann et al..,2002) . This 
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argument also negates the solar migration scenario. Fi- 
nally, a suggestion that the presence of ^°Fe in mete- 
orites proves that the proto-solar nebula was affected by a 
nearby supernov a explosion was recently found less likely 
(|William£j|2009l) . 

(c) Early dust separation. One possible explanation is 
that the Sun was formed in an interstellar cloud where 
dust had been blown away beforehand, e.g. by lumi- 
nous massive stars. Another possibility could be that 
the late formation process of the Sun included accretion 
of gas that was metal-poor due to dust separation by ra- 
diation pressure, although it is questionable whether the 
accreting proto-Sun was efficient enough in this respect. 
Neither of these scenarios can explain why this happened 
to the Sun but not the twins nor can it explain the cor- 
relation with presence of close-in giant planets (§4e). 

(d) Dust separation during the formation of terrestrial 
planets: A fascinating possibility is that the composi- 
tion departures of the Sun are related to its properties 
as a planet host. A particularly striking circumstance 
is that the inner solar system planets and met eorites 
are enriched in refractories compared to volatiles (jPalmd 
[2000). The abundance pattern of these bodies is almost 
a mirror-i mage of the solar pattern relative to the so- 
lar twins (|Cieslal 120081 : [Alexander et al.ll2001f) . Also, a 
radial gradient exists, with greater enhanc ement of re- 
fractories a t smaller heliocentric distances ('P almd [20001 : 
IBottke et~ al. 2006). The break in the chemical abun- 
dance pattern at Tcond ~ 1200 K (Fig . 3), suggests 
that the volatiles retained their original abundances both 
in the Sun and the solar twins. Such temperatures are 
only encountere d in the inne r parts {< 3 AU) of proto- 
planetary disks (jCieslal [20081 ). which also suggests that 
the abundance pattern is related to the presence of ter- 
restrial planets. 

The amount of dust-depleted gas required to explain 
the solar abundances depends sensitively on the timing 
of the accretion. The mass fraction of the solar con- 
vection zone dropped from 100% to 38% in 10 Myr, 
and to 2% after another 20 Myr , there after remaining 
largely constant (jP'Antona k, M azzitcUi 1994). Assum- 
ing that the depletion signature was imprinted once the 
convection zone reached its present size, the removal of 
~ 2 X 10^^ g of refractories from the accreting gas would 
be required in order to accomplish the observed ^ 10% 
reduction of refractory elements. This value is similar to 
the combined mass of refractories locked up in Mercury, 
Venus, Earth and Mars, which is 8 x lO^^g (LSM^). 
The formation of large (~ 10 km) bodies and embryos 
of terrestrial planets in the prot o-solar disk is estimate d 
to have an efficiency of 30-50% ([Alexander et al.ll200lD . 
Thus, it is tempting to speculate that the formation of 
the terrestrial planets might have given the Sun its spe- 
cial surface composition. The disk masses during the 
T Tauri phase are '-^ 0.02 M p, but values up to IMq 
are possible ([BeckwithI [l990l: iCalvet et all I2OOO0 . thus 
enough material would be available to change the solar 
photospheric composition. However there is a problem 
with time-scales since proto-planetar y disks are observed 
to have typical life-times < 10 Myr (jCalvet et all 120001 : 
ISicilia-Aguilar et "al1l2006[ ). Yet, gas accretion onto 10- 
25 Myr old stars has been detected but is rare: only 
1% of stars at an age of 13 Myr sh ow signs of accre- 
tion ([Currie et al][2007l: iWhite fc Hillebrand.2005.1 . One 



possibility to explain the chemical differences with the 
solar twins is that they did not form terrestrial plan- 
ets or, if they did, the subsequent gas accretion ocurred 
when the stellar convection zone was still deep and the 
planetary signatures thus washed out, while the Sun re- 
tained its planetary disk longer. A less likely explanation 
is that all solar twins formed terrestrial planets and like 
the Sun lost their gas disks when their convection zones 
were deep, but accreted their terrestrial planets later (af- 
ter 30 Myr) when the convection zones were thin, thus 
enriching the solar twins in refractory elements. There 
are reasons to believe that thermohaline convection, gen- 
erated by the inverse yit-gradient, would then have quickly 
erase d the higher m etallicity of the stellar convection 
zone (|Vauclaiill200l . 

(e) The role of giant planets: The abundance pattern 
could also be the result of the formation of giant plan- 
ets. The cores of Jupiter, Saturn, Uranus and Neptune 
containing heavy metals in the form of rocks and ices are 
estimated to have a combined mass of - 30Me (IGuillotI 
120051 ). Thus the accretion could have occurred earlier, 
which would alleviate the problem with disk life-times. 
If so, stars with giant planets would be expected to show 
a chemical abundance pattern more similar to solar. Our 
sample of 10 solar analogs from planet surveys includes 
four solar analogs with giant planets, while for the other 
six no planets have been detected yet. As shown in Fig. 
5 all solar analogs with giant planets differ from the Sun 
but closely resemble most solar twins. Thus, the odd so- 
lar composition is not due to giant planets as such. On 
the other hand, only two of the six solar analogs without 
close-in giant planets have abundances that differ signif- 
icantly from the solar pattern. The fraction of stars with 
the solar pattern seems thus tentatively related to the 
presence of giant planets on close orbits: ~ 0% when 
having such planets, « 20% for solar-type stars in gen- 
eral and w 50 — 70% without close-in giant planets. No 
doubt, the statistics has to be improved considerably, but 
the numbers are clearly tantalizing. 

5. CONCLUSIONS 

We conclude that Sun is unusual albeit not unique in 
its detailed elemental abundance pattern. In particular 
it has a chemical composition that is more affected by 
dust condensation than most other similar stars, such 
as most solar twins and all solar analogs with discov- 
ered giant planets. This may be purely circumstantial. 
For instance, the cloud in which the Sun was formed 
may first have been cleansed from dust by radiation 
from hot luminous stars to a higher degree than hap- 
pened where most other solar-type stars were formed. 
However, if an explanation of less ad-hoc character is 
searched for, it is natural to explore whether the solar 
composition could be related to the way the formation 
of the planets in the solar system occurred. A reason- 
able hypothesis is then that the accretion onto the Sun 
of the proto-planetary solar nebula, chemically affected 
by dust condensation and planet formation, was delayed 
long enough for the cleansed gas to significantly affect the 
composition of the solar convection zone. Another pos- 
sibility is that the early Sun was never fully convective, 
which would make it easier to imprin t a dust-cleansed 
abundance signature ([Nordlundl |2009[ ). This scenario 
is supported by the dynamical star-formation calcula- 



4 



Melendez et al. 



tions with a time-dependent conv ection treatment of 
iWuchterl fc Tscharnute^ (|2003f) and lWuchterl fc KlessenI 
(|200in . 

The fact that stars with known close giant planets do 
not show the solar pattern suggests that their gas disks 
for some reason were accreted earlier when the stellar 
convection zones were so deep that the planetary signa- 
tures in the infalling gas disks were erased, and/or that 
the key factor behind the solar pattern is the formation 
and non-accretion of inner terrestrial planets. This would 
imply that solar-like stars with planetary systems simi- 
lar to our own are a relatively rare occurrence. Before 
firm conclusions can be drawn regarding the effects and 
frequency of planet formation, accurate analyses of sig- 
nificantly larger samples of solar-like stars with and with- 
out known giant planets are important. It is reassuring 
that an independent analysis of 22 northern solar twins 
confirm our findings both in terms of the magnitude of 
the abundance difference between volatiles and refracto- 
ries and the frequen cy of stars showing the peculiar solar 
abundance pattern (| Ramirez "eFaLl [20091. Our findings 
provide strong impetus to improve the abundance accu- 
racy of late-type stars to the 0.01 dex level by remov- 
ing the remaining systematic errors in t he modelling of 
stellar atmospheres and line formation (jAsplundl 120051 : 
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Fig. 1.— High quality {S/N ~ 450) MIKE spectra of the Sun and solar twins in the region 5970-6030 A. 
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Fig. 2. — Differences between [X/Fe] of tlie Sun and the mean values in the solar twins as a function of atomic number Z. For clarity 
the elements Y (Z = 39), Zr (Z = 40) and Ba (Z = 56) have been included after Zn. Observational la errors in the relative abundances 
(including observational errors in both the Sun and solar twins) are shown with dotted error bars, while the Icr errors in the mean abundance 
of the solar twins are shown with solid error bars. 
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Fig. 3. — Differences between [X/Fe] of the Sun and the mean values in the solar twins as a function of T^ond- The abundance pattern 
shows a break at T^ond ~ 1200 K. The solid lines are fits to the abundance pattern, while the dashed lines represent the standard deviation 
from the fits. The low element-to-clement scatter from the fits for the refractory (tr = 0.007 dex) and volatile (cr = 0.011 dex) elements 
confirms the high precision of our work. The zero-point for the differences in relative chemical abundances depends on the adopted reference 
element, which here is Fe; the volatilcs would appear normal while the refractories more depleted had we instead selected to use C. Error 
bars as in Fig. 2. 
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Fig. 4. — Mean [X/Fe] ratios observed in the solar twins for highly volatile (triangles), semi-volatile (squares), medium refractory 
(pentagons) and highly refractory (stars) elements, for the solar twins. Most twins show similar trends with Tcond with the highly 
refractory elements (1500 < T^ond < 1750 K) well separated from the highly volatile elements {T^ond < 500 K) and in between the 
semi-volatile (500 < T^ond < 1250 K) and medium refractory (1250 < T^ond < 1500 K) elements. This separation between refractories and 
voiatiles on a star-by-star basis further emphasizes that the abundance differences are indeed real. 
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Fig. 5. — Differences between [X/Fe] of the Sun and the mean values in the solar analogs with (filled triangles) and without (open squares) 
detected planets, as a function of their Tcond- The planet hosting stars show a strong trend with Tcond almost identical to that seen in 
the solar twins (Fig. 3). The correlation is highly statistically significant (rg = —0.91) with a probability of having this correlation by 
pure chance of ~ 10~*. Solar analogs without planets show a much shallower slope and a much weaker correlation. Linear fits to both 
samples are shown. The slope of the stars with planets is significant at the 7cr level, while the slope of the stars without planets is lower 
and significant only at the 2ct level. 



